Quantum Spherical Spin Glass. 
Supersymmetry and Annealing 

Pedro Castro Menezes 
Alba Theumann 
Institute de Ffsica, Universidade Federal do Rio Grande do Sul 
Av. Bento Gongalves 9500, 91501-970 Porto Alegre, RS, Brasil 



Abstract 

We show that the effective action of the quantum spherical 
spin glass is invariant under a generalized form of the Becchi- 
Rouet-Stora-Tyutin (BRST) supersymmetry. The Ward identi- 
ties associated to this invariance indicate that the spin glass order 
parameter must vanish, and as a result the annealed average is 
exact in this model. We present new results for the free energy, 
entropy and specific heat. Due to quantum effects the entropy 
remains finite and the specific heat vanishes at zero temperature. 
Results for the phase diagram coincide with those obtained by 
different formalisms. At zero temperature we derive the scaling 
behaviour with frequency of the dynamical susceptibility . 



PACS numbers: 75.10.Nr; 75.10.Jm ; 64.70.Pf 
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1. Introduction 



In a previous publication [I] the classical spherical spin glass was analyzed by 
using super symmetry methods [21 E] and it was pointed out that the model 
was "anomalous" because the order parameter q that measures the replica's 
overlap was vanishing and the configurational average of any number of repli- 
cas gave the correct known result [1], then the annealed approximation with 
only one replica was also correct. More recent work on the calculation of 
the complexity, or logarithm of the density of states of a certain energy, of 
the p-spherical spin glass, p > 3, [5j [6] and of the Sherrington-Kirkpatrick 
model [7] shows that the annealed average gives the correct answer when the 
action is invariant under the generalized Becchi-Rouet-Stora-Tyutin (BRST) 
supersymmetry[8] . 

In the present paper we first demonstrate how the generalized BRST super- 
symmetry transformation applies to the functional integral formulation of 
the quantum spherical spin-glass. We also discuss how this internal invariance 
leads to Ward identities that result in the vanishing of the replica overlap 
g.This result had been obtained previously for the classical spin-glass[T] and 
justifies the performance of the annealed configurational average [51 El Ej also 
in the quantum spherical spin-glass. The model was proposed earlier in ref. [9] 
and was solved by using the distribution of eigenvalues of a gaussian random 
matrix. Later, the interest in this model was renewed because it corresponds 
to the infinite-M limit of the spin-glass of M-components quantum rotors 
and it was studied [10] by using the replica method in a replica symmetric 
theory. Both theories give the same phase diagram, although thermody- 
namic quantities as free energy and specific heat are not discussed. On the 
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other hand, the spin-glass model of quantum rotors is equivalent to the Ising 
spin-glass model in a transverse fieldJTT]. A formulation for quantum spher- 
ical spins with a different temporal dependence was presented for multispin 
interactions (12]. Critical properties were analyzed in ref. [13] and applications 
to physical systems were presented in ref. [TH [To] . All these papers coincide 
in that the quantum spherical model presents a finite entropy and vanishing 
specific heat at low temperature, while the classical model presents a nega- 
tive infinite entropy and constant specific heat. 

We study the quantum spherical spin-glass in the annealed average by us- 
ing a functional integral formulation derived by using Feynman's prescription jT6j 
ITT] . In sect. 2 we present the general formalism and discuss the supersym- 
metry transformation that leaves invariant the action together with the cor- 
responding Ward identities [18J, while in sect. 3 we discuss the results. We are 
careful to normalize the free energy so as to recover a finite classical limit 
when the rotors moment of inertia I — > oo. In this way the phase dia- 
gram, that is calculated from the spherical condition, coincide with previous 
studies jH [10]. In the classical limit / — > oo the entropy diverges at T = in 
the typical fashion of the classical spherical model[TJ H], with the important 
difference that it remains finite in the quantum regime for finite values of the 
rotors moment of inertia. The specific heat shows the usual discontinuity at 
the critical temperature but it vanishes at T = in the quantum regime. We 
derive scaling laws for the dynamical susceptibility, in agreement with pre- 
vious results [10]. We present in the Appendix some useful relations among 
Grassmann variables without proof, while we refer the reader to the original 
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work[2] for rigorous results. 



2. Supersymmetry and Annealing 

We consider a spin glass of quantum rotors [10] with moment of inertia I in 
the spherical limit with Hamiltonian 

j i i,j i 

where the spin variables at each site are continuous — oo < Si < oo and we 
introduced the canonical momentum P, with commutation rules: 



[Sj,Pk]=i6M (2) 

The sum in eq.((T])runs over sites i,j = 1..JV. The coupling in eq.flT]) is an 
independent random variable with the gaussian distribution 



P(J«)=e-^^ (3) 

while the chemical potential /x is a Lagrange multiplier that insures the mean 
spherical condition 

dilnZ) „ /•/* 



W dT($) = PN (4) 
„■ jo 



and /3 = 1/T is the inverse temperature. We work in units where the Boltz- 
mann constant kn = h = 1 and Q is the quantum partition function 

z = Tre -m SG +^^) ( 5 ) 



The partition function at zero field can be expressed as a functional integral [TBI 

Q3] 



z 




(6) 



where the non interacting action Ao is given by 



Ao 
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+ ^St{r)-H i {r)S i {r) 



) 



(7) 



and the interacting part 



Asg = I] Jij / dTSi(T)Sf(r) 




(8) 



After performing a Fourier transformation in time and taking into account 
the reality of the fields Si(r) the generating functional can be written in terms 
of the real components Ri(n) of the Fourier transforms Si(uj n ) = S*(—u n ), 
where u n = are boson Matsubara's frequencies. As it was shown in [1] 
the correlation functions can be obtained from the generating functional for 
two replicas 



w(#) = nn r *;(™)ii /.n dRiain) 



n>0 ian>0 



exp - Y^ij{n)R ia {n)R ja (n) + Y,H ia (n)R ia (n) 




(9) 



where the index a = 1, 2, the Hi a (n) are auxiliary fields and 



T i3 (n) = {1(3^1 + 2(3^-3, 



(10) 
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By using the results in the Appendix , the determinant ||rVj(n)|| may be 
expressed with the help of auxiliary Grassmann fields x*^X an d W can be 
written 

W(H a ,-y) = II / II d X :(n)d Xt (n)Y[dR ia (n) 

n>0 in>0 a 

X*(n)xj(n) + X*j(n)Xi(n) 



exp 



n>0 ij a 

+ 2^ o + 2^H iQ {n)R ia {n) 



where Xi( n )tXi( n ) are complex anticommuting Grassmann variables while 
the Ri >a (n),a = 1,2 are real commuting variable and we introduced two 
extra Grassmann auxiliary fields 7$ and 7*. As the T^- are symmetric, when 
the auxiliary fields are set equal to zero the functional is invariant under the 
supersymmetry transformation 



x£»/V2 = X-W/V2 + e*XX a (n)/V2 
Xi(n)/V2 = Xi(n) / V2 + eJ2 R iA n ) / ^ 



ifi, a (nj = i£i^(n) 1 — (12) 

where e is a complex Grassmann variable and we adopt the convention for 
complex conjugation [2] e** = -e, (ex*)* = e*X** = ~ e *Xh then R ita (n)' 
is effectively a real variable. The transformation differs infinitesimally from 
unity and e*e can be neglected, although this concept is somehow meaningless 
in the case of Grassmann variables. Now it is convenient to switch to the 
superalgebra notation [2 J of the Appendix and to introduce the supervectors 
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for the auxiliary fields 



SJ] 



( H a \ 

H a H l2 7*/V2 -7i/V2) 



(13) 



and for the field variables 

Xi/V2 

\xVV2J 

S^ t = (Ra Ri 2 X1/V2 -Xi/V2) (14) 

In this way the functional W in equation ffTTl) may be written in compact 
form 



W(J) = II / II dx*(n)d Xt (n)l[dR ta (r 

^1 in>0 a 



cxp 



1.1 



that may be formally integrated to give 



W(J) = exp 
with the 4x4 supermatrices 



1/2 X^Qii Jj 



(15) 



(16) 



Qij 



/ Qll <?12 #13 #14 \ 

921 <?22 #23 #24 

#31 #32 ^33 

V#41 #42 VteJij 
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(17) 



By symmetry we should have 



qil = <?22 =< RiaRia >= V 33 = V U =< X* X > 

v 34 = v 43 =< XX >=< X*X* >= 



(18) 



To derive Ward identities [18] we consider that, as the system is invariant 
when the fields are subject to the supersymmetry transformation in equa- 
tion^ 



where 



/ 1 





1 



1 + D 



(19) 



-e*/V2 e/ v / 2\ 
-e*/V2 e/y/2 
e/V2 e/y/2 1 
\e*/V2 e*/V2 1 / 

then it should be also invariant when we transform the external fields [18] in 
equation tfTfi]) 



J' = A J 
SJ P = SJ+SA* 



(20) 



Where the super-adjoint SA^ is defined in the Appendix. This condition 
leads to 



SD f Q + QD = 



(21) 



from where we deduce 



^33 = qn + qu 

ViA = <?22 + <?21 
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(22) 



and from equation ([TS]) we obtain for the overlap between two replicas 

<?i2 = g-si = (23) 

The spin-glass order parameter q = g 12 and it vanishes identically from equa- 
tion f|T9|) . thus making the annealed average exact in the quantum spherical 
spin-glass as it was discussed in[T] for the classical spherical spin-glass. 



3. Results 



When performing the annealed average we take the configurational average 
of Z in equation(j6]) over the random variables Jy and after splitting the 
quadratic term with a gaussian integration we obtain the result 

J 2 f p 



(Z)ca = J ~ r') exp -N J q jf drdr'Q{r - r 

where we assumed time translational invariance and 

J 2 fP rP 



f\2 



A 



(24) 



VS{t) exp 



-*4o + V / / drdr'Q(T - t')S(t)S(t') 
l Jo Jo 



(25) 



A steepest descent calculation in eq. pi)) gives 

Q(r-r') = ^Y:<Si(r)S i (r f )> 



(26) 



The partition function in eq. (!24p is solved by a time Fourier transformation 
with the result 

■(/3J)*, - /3H*(u n )H(u n ) 



{Z)ca = ^ II / d Q^n) exp{-iV[^QK) 2 - 

n=0 J 1 



(3Ioo n 2 + 2/3/i - 2{(3J) 2 R{uj n ) 
+ ln (piu 2 n + 2Pfi-(pj) 2 Q(u; n ))]} 



(27) 



where uo n = is the boson Matsubara frequency and Q(u n ) is the real 
part of the Fourier transform of Q(r). A steepest descent evaluation of the 
integral in eq. fl271) gives, for H(ui n ) = 

Q(u n ) = [f3Iu 2 n + 2$n - {l3J) 2 Q{u n )}- 1 (28) 

with solution 



/3J^ + 2 /3 / ,- v /(^ + 2/? / i)2-4(/5J) 2 
QM = ^ (29) 

while the mean spherical condition in eq. ([l"l) becomes 

oo 

£QK) = i (30) 



Introducing eq. ff29|) into eq. ff27j) we obtain for the free energy per site at the 
saddle point 

P F = £ (j{PJ?Q{un) 2 - ^ In (QM) J - - In (Af) (31) 
The normalization constant M = —3= is determined so that we recover a finite 



limit for the free energy when / — > oo . Following standard procedures [16] 
we convert sums over frequencies into integrals with the result for the mean 
spherical condition 

f L + dyy/(Ll-y*)(y*-L±)coth(^L) = 2ttJ 2 V7 (32) 

where 

L 2 ± =2fi± 2 J (33) 
and we obtain for the free energy 

PF = VP Il + rf ^ lnsinh V / 4J 2 -(2/i-y 2 ) 2 -fin- \nM 

(34) 



In the classical limit / — > oo the integrals in eq. (1321) and eq. (1341) can be 
performed exactly in terms of hypergeometric functions with the result 

^f^ = pj+jj (35) 



I , . — — 1 fa - v/(/3/i) 2 - {(3 J? 

+ T , = - HP 

4/?/x+^/(/?/x) 2 -(/?J) 2 



= - In [ >'// + v/(^) 2 - (/3 J) 



2 



2 

(36) 



The expression in eq. f[35|) coincide with the results for the classical spherical 
spin glass in ref . ( [Tj H] ) , but the free energy in eq. (l3"61) differs due to the last 
logarithmic contribution. This originates in the choice of the normalization 
constant Af. If we had chosen M following Feynman's prescription |16[ [T7] we 
would have obtained unphysical results in the quantum regime for finite val- 
ues of /, like a negative specific heat at low temperatures [T9] . The chemical 
potential fi decreases with temperature until it reaches the value /i c = J[20\ 
at the critical temperature T c , for lower temperatures the spherical condition 
no longer holds and fi sticks to this value. We show in fig.l the solution 
for fi as a function of T for different values of /. The critical temperature 
is obtained by setting \i = // c = J in eq. (!35|) . We show in fig. 2 the phase 
diagram for T c as a function of J. T c decreases with 1/J and vanishes at the 
critical value I c = g^fj. This result coincides with previous references (9j ITU] . 
Close to this value we find T 2 = ^ "^j/^ • We show in fig. 3, fig.4, fig. 5 
results for the free energy, entropy and specific heat as functions of temper- 
ature for different values of I. We observe that the entropy remains finite 
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when T — > for any finite value of J, while we recover the usual divergence 
S — > — oo in the classical spherical spin glass modelp], H] when I — > oo. The 
specific heat vanishes at T = for finite values of /, as expected in quantum 
systems [T2l [TBI ITU [15]. The dynamic susceptibility is the response to the 
field H{u n ) and we obtain from eq. (1271) 



^/ + 2/i- v /(^/ + 2 /U ) 2 -4J2 
XK) = ^775 (37) 



2J 2 

The static susceptibility 



*(°) = ( 38 ) 

satisfies Curie law for high temperatures and exhibits the usual cusp at the 
critical temperature. It is shown in fig.6 for different values of /. The dynamic 
susceptibility at T = is obtained by setting iuj n — >• uj + i8 in eq. ( 1371) .For 
A = jj, — J > and I = I c we obtain that the Imx{oj) satisfies the scaling 
relation 



Im X (v) = sgnulul'iy/TjJ)*^ -j^) (39) 

with the exponent // = 1 and the scaling function $(x) = yl — a; -2 , what 
agrees with the results of ref.([XU]). 



4. Conclusions 

The quantum spherical spin glass was studied within the annealed average, 
that is exact due to the internal BRST supersymmetry of the model, as it 
was discussed in sect.l. Our results for the phase diagram coincide with 
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those obtained by Ye et alfTU] with the replica method and by Shukla et 
al [9] using the semicircular law. However, a substancial difference exists 
between our results for the free energy, entropy and specific heat and those 
presented by Shukla et al[9]. The results obtained with quantum functional 
integrals methods are very sensitive to the use of the correct normalization 
and measure [T6l [T7] and we set the normalization constant in our expression 
for the free energy such as to recover a finite limit when I — > oo. As a 
result we obtain a finite entropy at low temperatures and a vanishing specific 
heat for finite values of the rotors moment of inertia /, as it is expected in 
quantum spherical models [121 ESI US]- I n the limit / — ► oo we recover the 
known negative, infinite entropy and a constant specific heat at T = 0. Also 
at T = we derive the scaling behaviour for the dynamical susceptibility 
that agrees with Ye et alfTO]. 
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5. Appendix 



We refer here to references [2J [3] to give a concise and brief description of 
Grassmann variables and superalgebra. We indicate complex Grassmann 
variables by Greek letters, and these anticommuting variables satisfy 



[x,7]+ = X 2 = X 



* 2 







x = -x 
(xyY = x*i* 

\X X) = X x = -xx = X X 



(40) 



and the integrals 



what leads to 



Jd X = Jd X * = 
J U d XidXiexpJ^xtM ijXj = \\M K 



(41) 



(42) 



A four-component supervector <p and a 4x4 super matrix M have commuting 
and Grassmann components 

/ Ri \ 

R2 

X/V2 

W/V2J 
a a 

P h 



M 



(43) 
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where a and b are 2x2 matrices of commuting elements while g_ and p are 
2x2 matrices of Grassmann elements. Taking into account the change of 
sign involved in complex conjugation from eq (j4"0]) the super adjoint must be 
defined 

V = (^ r 2 r /^2 -x/y/2) 

5M f = "I (44) 



14 



6. Figure Captions 



Fig.l Chemical potential /i as a function of temperature for different val- 
ues of the moment of inertia / = 0.19J(dash-dot), I = J (dash), / = 
oo (continuous) . 

Fig. 2 Phase diagram and critical line T c (l/J) separating the paramag- 
netic from the spin glass phase in the Tvsl/I plane. 

Fig. 3 Free energy as a function of temperature for different values of 
the moment of inertia / = 0.4J(dash), I = J (dash- dot- dot), J = 5J(dash- 
dot),/ = lOJ(dash-dot-dot) I = oo (continuous) . 

Fig. 4 Entropy as a function of temperature for different values of the 
moment of inertia / = J (dash), I = 5J(dash — dot — dot), I = 10 J (dash — 
dot), I = oo (continuous) . 

Fig. 5 Specific heat as a function of temperature for different values of the 
moment of inertia / = 0.4 J (dash),/ = J (dash-dot), / = 10J(dash-dot- 
dot),/ = oo (continuous) . 

Fig. 6 Static susceptibility as a function of temperature for I = 0.4 J(dash), 
I = J (dash-dot),/ = oo (continuous) . 
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